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Abstract. We consider the static Holstein model, describing a chain of 
Fermions interacting with a classical phonon field, when the interaction 
is weak and the density is a rational number. We show that the energy of 
the system, as a function of the phonon field, has two stationary points, 
defined up to a lattice translation, which are local minima in the space 
of fields periodic with period equal to the inverse of the density. 
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1. Introduction 

1.1. The Holstein model [P,H] was introduced to represent the interaction of electrons with 
optical phonons in a crystal. In the original model the phonons are represented in terms of 
quantum oscillators but the difficulty of understanding such a fully quantum model has led 
to a modification of it, called static Holstein model (or adiabatic Holstein model), in which 
the phonons are classical oscillators. This corresponds to neglect the vibrational kinetic 
energy of the phonons, an approximation which can be justified in physical models as the 
atom mass is much larger than the electron mass. 

The Hamiltonian of the model, if we neglect all internal degrees of freedom (the spin, for 
example, which play no role at zero external magnetic field) is given by 

xeA 

x.yEA i£A xEA xEA 

where x, y are points on the one-dimensional lattice A with unit spacing, length L and 
periodic boundary conditions; we shall identify A with {ieZ: — [L/2] < x < {{L — l)/2]}. 
Moreover the matrix t xy is defined as t xy — S XtV — (l/2)[<5 Xi!/+ i + <5 Xi!/ _i], where 5 XiV is the 
Kronecker delta, p is the chemical potential and A is the interaction strength. The fields 
ip x are creation (+) and annihilation (— ) fermionic fields, satisfying periodic boundary 
conditions: Vx = i'x+L- Wc dcfine also V>x = e tH i))±e~ Ht , with x = (x, t), -0/2 < t < 0/2 
for some > 0; on t antiperiodic boundary conditions are imposed. The potential <p x is a 
real function representing the classical phonon field. 

At finite L, the fermionic Fock space is finite dimensional, hence there is a minimum 
eigenvalue Ef(<p, p) of the operator H L l , for each given phonon field ip and each value of p; 
let pl((P,[J.) be the corresponding fermionic density. The aim is to minimize the functional 

F L (<p,n)=Et(cp,n) + ±J2vl> (i- 2 ) 

subject to the condition 

Pl(<P,») = PL , (1.3) 

where pl is a fixed value of the density, converging for L — > oo, say to p. 

The model (1.1) can be considered as an approximation of a more realistic continuous 
model containing also the interaction with a fixed external periodic potential of period one. 
Then the discreteness is not a pure mathematical artifice, but it has a precise physical 
interpretation: the properties of the two models are expected to be the same, and we think 
that this could be easily proven along the lines of the present paper. 

1.2. It is generally believed that, as a consequence of Peierls instability argument, [P,F], 
there is a field <p(°\ uniquely defined up to a spatial translation, which minimizes (1.2), 
(1.3), in the limit L — > oo, and is a function of the form <p(2npx), where (p{u) is a 2-k- 
periodic function. This is physically interpreted by saying that one-dimensional metals are 
unstable at low temperature, in the sense that they can lower their energy through a periodic 
distortion of the "physical lattice" with period 1/p (in the continuous version of the model, 
since 1/p is not an integer in general): such a distortion is called a charge density wave, 
as the physical lattice and the electronic charge density form a new periodic structure with 
period bigger than the original lattice period. 

The argument in [P,F] is quite simple: the periodic potential (p{2-Kpx) opens a gap in 
the electronic dispersion relation in correspondence of the Fermi momentum, and a trivial 
computation using degenerate perturbation theory shows that the elastic energy increase 
is less than the fermionic energy decrease. However, see [LRA], in this argument one does 
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(1.1) 



not take into account the effects due to the discreteness of the lattice, in particular the 
fact that the momenta are conserved modulo 2tt ( Umklapp) . Neglecting the discreteness of 
the lattice one loses the difference between commensurate or incommensurate charge density 
wave (i.e. rational or irrational p) in the infinite volume limit, whose properties are supposed 
to be different, especially concerning the conductivity [F,LRA]. 

Note also that, even if the argument in [P,F] is perturbative, Peierls instability is expected 
to arise also for large interaction strength, [AL] . 

An exact result, [KL,LM], makes rigorous the theory of Peierls instability for the model 
(1.1) in the case p = pl = 1/2 (half filled band case), for any value of A. In fact, in this case 
it has been proved that there is a global minimum of Fl(<p) of the form e(X)(—l) x , where 
e(X) is a suitable function of A. This means that the periodicity of the ground state phonon 
field is 2 (recall that in our units 1 is just the lattice spacing): this phenomenon is called 
dimcrization. The proof heavily relies on symmetry properties which hold only in the half 
filled case. 

In [AAR,BM] Peierls instability for the Holstein model is proven assuming A large enough: 
in that case the Fermions are almost classical particles and the quantum effects are treated as 
perturbations. The results hold for the commensurate or incommensurate case; in particular 
in the incommensurate case the function <p(u), related to the minimizing field through the 
relation ip x = {p(2iTpx), has infinite many discontinuities. On the contrary, in the small A 
case, according to numerical results, (p(u) has been conjectured to be an analytic function 
of its argument, both for the commensurate and incommensurate cases, [AAR]. 

In this paper we study the case of small A and any rational density, for which there are, 
to our knowledge, no results in the literature besides the simulations in [AAR]. Analytical 
results in the small A case can be found for a related model, the Falikov-Kimball model, 
described by a Hamiltonian of the form (1.1), in which the continuous (p x is replaced by a 
discrete function taking only the values or 1; see [FGM]. 

Let p = P/Q, with P, Q relative prime integers, and let L = Li = iQ, i € N; we shall 
prove that, if A is small enough, there are two stationary points tp( ± ' % ') of Fz^ (<p, /x), defined 
up to a lattice translation, satisfying (1.3) with pl — p. These stationary points are periodic 
functions on L{ of period Q (the smallest multiple of 1/ p which is an integer, hence a multiple 
of the unit lattice spacing), converging for i — > oo. Moreover, if we restrict Fl^^p) to 
functions such that (p x = <p x +Q, (fi^'^ are local minima in the norm \\tp\\ = sup.j, ei |<p x |, 
uniformly in i. 

The presence of the lattice has the effect that we need the smaller A the bigger Q is, see 
(1.16). In particular we are not able to draw conclusions about the incommensurate case 
neither we know if this is a technical limitation or there is some physical reason behind it, 
so that we can not draw any conclusions about the analyticity conjecture in [AAR]. 

1.3. Let h xy — t xy — \(fi x 5 xy be the one-particle Hamiltonian and ei((p) < e2(<p) < ■ ■ ■ < 
e z,(<£) its eigenvalues. We have 



where P M is the projector on the subspace spanned by the eigenvectors of h with eigenvalues 
< p. As it is well known (see, for example [BM]), Ef((p 7 p) is a differentiable function of <p 
and, since TiChdP^/dipx) = 0, 



where p x (ip,p) = (P^)^ is the density of the electrons in the point x. 

Let us now suppose that p is not equal to any eigenvalue of h. In this case, given (p, 
also Pi((y5, p) is differentiable in a neighborhood small enough of <p (so small that e n (ip) — p 





(1.4) 



n:e n (ip)<n 




(1.5) 
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stays different from zero, for any n, see again [BM]) and dpL{p 7 p)/dp x = 0. Hence, a local 
minimum of (1.2) satisfying (1.3) must satisfy the conditions 

fx = Xp x (f,p) , 

PL = J^Px{V,IJ) , < ~ 1 ' 6 ' ) 



M xy = S xy — X— — p y (p 7 p) is positive definite . (1.7) 



d_ 

'df x 

Note that, given if, the previous condition on p can be in general satisfied only if | \<p — if\ \ 
is of order l/L, so that a solution of (1.6) defines in general a local minimum only in a 
neighborhood of size l/L. It follows that the only solutions which are interesting in the 
limit L — > oo are those associated with a gap of h around p, whose size is independent of L. 

1.4. If ip is a solution of (1.6), it must satisfy the condition ifo = ip x — XpL- On the 

other hand, if we define \x = f x —fo, we can see immediately that phif, p) = Pl(x, p+Xfo). 
It follows that we can restrict our search of local minima of (1.2) to fields ip with zero mean, 
satisfying the conditions 

fx = X(p x (f,p) - Pl) , 

PL = j^2px(<P,p) , ^ 1 ' 8 ' ) 

X 

and condition (1.7). 

Of course, if the field ip x satisfies (1.8), the same is true for the translated field p x+ni for 
any integer n. On the other hand, one expects that the solutions of (1.8) are even with 
respect to some point of A; hence we can eliminate the trivial source of non-uniqueness 
described above by imposing the further condition ip x = tp- x . We shall then consider only 
fields of the form 



[(L-l)/2] 

<Px= £'n e ^ . 0-n = f'n G K , ^0 = 0. (1.9) 

n=-[L/2] 



As we said in §1.2, we want to consider the case of rational density, p — P/Q, P and Q 
relatively prime, and we want to look for solutions such that ip x = f x +Q- Hence, we shall 
look for solutions of (1.8) with L = Li = iQ, pl = P and 

I(Q-i)/2] 

fx= Yl Vn^ pnx , p n = p- n eR, p = 0. (1.10) 

n=-[Q/2] 

Note that the condition on L allows to rewrite in a trivial way the field ip x of (1.10) in the 
general form (1.9), by putting (p' n = for all n such that (2nn)/L ^ 2i/)m,Vm, and by 
relabeling the other Fourier coefficients. 

The conditions (1.8) can be easily expressed in terms of the variables <f n ; if we define p n 
so that 

KQ-i)/2] 

p x {<p, f i)= p-n{f,py 2n * px , (l.n) 

n=-[Q/2] 

we get 

<p n = \ Pn {y, l i) , n^0, n = -[Q/2],...,[(Q-l)/2], (1.12) 

p (<f,fi)=p L . (1.13) 

Also the minimum condition (1.7) can be expressed in terms of the Fourier coefficients; 
we get that the L x L matrix 

d 

M nm = S nm - X— — p' m {f,p) (1.14) 
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has to be positive definite, if the field ip satisfies (1.12) and (1.13) and p' m (ip,p) is defined 
analogously to (p' m in (1.9). Hence, if we restrict the space of phonon fields to those of the 
form (1.10), we have to show that the Q x Q matrix 

d 

M nm = S nm - A— —p m (ip,p) (1-15) 
o<p n 

has to be positive definite, if the field ip satisfies (1.12) and (1.13). 

1.5. Remark. It is easy to show (by using the expansion described in §3, for example) 
that Mnm can be different from zero only if 2ir(n — m)/L is of the form 2wpk for some k. 
However, we are not able to get good bounds on all matrix elements M nm not belonging 
(up to a relabeling of indices) also to the matrix (1.15); therefore, in studying the minimum 
condition, we restrict ourselves to the fields of the form (1.10). 

1.6. Theorem. Let p = P/Q, with P, Q relative prime integers, L = Li = iQ. Then, for 
any positive integer N, there exist positive constants e, e, c and K, independent of i, p and 
N , such that, if 

- \og{ev L) " K"N\ \og(cQ/v%) ' 
where 

v = sin(Trp) , (1.17) 

there exist two solutions ip^ of (1.8), with L = L - t , l — p = cos(np) and pL = p, of the form 
(1.10). The matrices M corresponding to these solutions, defined as in (1.15), are positive 
definite. 

Moreover, the Fourier coefficients ip^ verify, for \n\ > 1, the bound 

Finally, Xip^ is of the form 

\-(±) +2 / 27 TO0 + /3W(A,L) 1 

\ip\ = ±v cxp I — j , (1.19) 

with 

\P^(X,L)\<CX 2 (l + log^) . (1-20) 

where C is a suitable constant. 

The one-particle Hamiltonian h corresponding to this solution has a gap of order \\<pi\ 
around p, uniformly on i. 

1.7. The above theorem proves that there are two stationary points of the ground state 
energy in correspondence of a periodic function with period equal to the inverse of the 
density, if the coupling is small enough and the density is rational, and that these stationary 
points are local minima at least in the space of periodic functions with that period. The 
energies associated to such minima are different so that the ground state energy is not 
degenerate. 

The theorem is proved by writing p x (ip, p) as an expansion convergent for small A and 
solving the set of equations (1.12) by a contraction method. As a byproduct we prove that 
the (p n are fast decaying, (see (1.18)), so that p x is really well approximated by its first 
harmonics (this remark is important as the number of harmonics could be very large) . 

The results are uniform in the volume, so they are interesting from a physical point of 
view (a solution defined only for |A| < 0(1/ L) should be outside any reasonable physical 
value for A). The case in [KL] for the half filled case is contained in Theorem 1.6, but in 
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[KL] it is also proved that the solution is a global minimum. On the other hand this case is 
quite special (see Remarks 2.5 in §2). 

Finally the lower bound in (1.16) is a large volume condition: this is not a technical 
condition as, if the number of Fermions is odd, there is Peierls instability only for L large 
enough. The upper bound for A in (1.16) requires A to decrease as Q increases: in particular 
irrational density are forbidden. This requirement is due to the discreteness of the lattice 
and to Umklapp phenomena. Note that the dependence of the maximum A allowed on Q is 
not very strong as it is a logarithmic one. 

The case of irrational densities (possible only in the infinite volume limit), excluded by 
our theorem, is physically interesting, but the existence of Peierls instability in this case is 
proven only for large A, [AL,BM]. In [BGM] p x (ip,fi) is shown to be well defined for small 
A not only in the rational density case, (in which the proof is almost trivial), but also in 
the irrational case: in fact the small divisor problem due to the irrationality of the density 
can be controlled thanks to a Diophantine condition. However to solve the set of equations 
(1.12) we use a contraction method which is not trivially adaptable in the latter case (see 
Remarks 2.5 in §2). The same kind of problem arises in proving the positive definiteness of 
M nm in the rational case (and this is the reason why we are able to prove that the stationary 
points are local minima only in the space of periodic functions with prefixed period). As 
we said above, we do not know if such problems are only technical or there is some physical 
reason for this to happen. 
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2. Solution of the self-consistence equation 



2.1. Let p = P/Q, with P, Q relatively prime integers such that < P < Q, and L = Li = 
iQ; we have to look for a solution of (1.12) and (1.13), which is well defined for |A| < sq, with 
£o independent of i (otherwise our solution is meaningless from a physical point of view). As 
discussed in §1.3, this means that our solution has to be looked for in a class of functions for 
which the one-particle Hamiltonian h has a gap around p of width independent of L. This 
class of functions is described by the following lemma, to be proved in §5.4. 

2.2. Lemma. Let ip x be a field of the form (1.10), L — Li, 1 — p, = cos(-Kp), \\(pi\ > 
and \X(p n \ < a\\ipi \/\n\ N for some positive constants a and N. Then there exists Eq > 0, 
independent of i and p, such that, if \Xipi\ < SqVq/Q, with v = sin(7rp), the one-particle 
Hamiltonian h has a gap of width > \Xipi\/2 around p. Moreover, p n (ip, p,) is a continuous 
function of X, which converges to a continuous function of X as i — > oo, and p (<p,p) = p. 

2.3. We can write the self-consistence equation (1.12) as 

(p n = -X 2 c n (cr)<p n + Xp n (a, $) , a = X(px, $ = {A<£„}|„|>i , (2.1) 

where c n (a) depends on ip only through a. We write p n as a perturbative expansion in A 
(different from the power expansion in A); this expansion is described in §3. If \n\ > 1, 
we are here defining —Xc n (a)(p n the contribution to p n proportional to ip n of order 1 in 
the expansion, while — aci(a) is the contribution to p\ proportional to a of order < 1 in 
the expansion (explicit expressions for c n (a) and Ci(a) will be given in (4.16) and (4.39) 
respectively); p n takes into account all the remaining terms of first order plus all terms of 
order higher than 1. 

The equation (2.1) has of course the trivial solution ip n = 0, Vn, but it is easy to see that 
this is not a local minimum, by using the expansion for p n of §3. Therefore we shall look 
for solutions such that a ^ 0, so that we can rewrite (2.1) as 

(1 + A^)) = ^^, (2.2) 
o~ 

(1 + X 2 c n (a)) 

Note that the equation for n = —1 does not appear simply because p_i = p\, as a conse- 
quence of the condition <p n = ip- n e R, see (1.10). 

2.4. We prove Theorem 1.6 in three steps as follows. 

• We first study the self-consistence equation (2.3), considering cr as a variable belonging to 
the interval 

J = ( -cxp(-7rv /A 2 ) , cxp(-7ru /A 2 ) ) • (2.4) 

We find a solution, that we denote $(<r), if A is small enough. 

• We then prove that, if L is large enough, the equation (in A) 

l + A 2 Cl (q) = ^^ (2 . 5) 
a 

has two solutions e J, of the form (1.19). Therefore (cr( ± '(A)/A, ^(a^{X))/ A) turn 
out to be solutions of (1.12), which verify, thanks to Lemma 2.2, (1.13) with L = Li. 

• We finally prove that the Hessian matrices (1.15) corresponding to these two solutions are 
positive definite. 

2.5. Remarks. The coefficient (p\ has a privileged role with respect to the other coefficients. 
In fact, as we shall see in §5, the properties of the system when only (p\ is different from 
are very close to the properties of the case in which all the coefficients are non vanishing. 
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This suggests that the "important" equation is (2.2), so explaining the strategy outlined 
above. 

The previous remark also implies that 1 + A 2 ci((t) ~ 0. It follows that 1 + X 2 c n (a) ~ 0, for 
all n such that 2irpn ~ 2np (mod27r). Since min^i-^ \2irpn — 2np\ = 2n/Q, we can expect 
that our bounds will not be uniform in Q. This is the reason why Theorem 1.6 can not be 
extended to irrational density; at most one can hope that a Diophantinc condition on p is 
needed, but we have only been able to prove that the Q dependence can be substituted with 
a dependence on the Diophantinc constants in some of the bounds described below. 

Note also that, if Q = 2, the only equation to discuss is just the equation (2.2) with $ = 
and the r.h.s. equal to zero; its solution is well known in this case, see [KL,LM] for example. 
If Q = 3, again (2.2) is the only equation to discuss, but the r.h.s. is different from zero; 
however it is easy to prove that the solution has essentially the same properties as in the 
case Q = 2. Hence, in the following we shall consider only the case Q > 4. The following 
lemma, furnishing a bound on the constants c„(cr) and their derivatives, is proven in §4.9. 

2.6. Lemma. There exists a constant C, independent of i and p, such that, if \n\ > 2, 

\c n (a)\ <-(l + log-) log Q , (2.6) 

%^ < ^ , (2-7) 
da v \a\ 



2.7. Fixed L = L i: $ is a finite sequence of Q — 3 elements, which can be thought as a 
which is a function of a. In order to study the equation (2.2) for cr, we 



vector in 



shall need a bound on $ and on the derivative of $ with respect to a. Hence we consider 



the space T = C 1 ( J , 1 



jQ-3 



) of C 1 -functions of cr e J with values in 



the solutions of 



(2.3) can be seen as fixed points of the operator Ta : T — > T, defined by the equation: 
We shall define, for each positive integer N, a norm in T in the following way: 



= sup <^ 
|n|>l,creJ I 



n 



iN 



\o-\- 1 \$ n (a)\ + 



da 



We shall also define 



B = {$ G T: < 1} ; 

R($) n (a) =p„(<7,$(<7)) , |n|>2. 



(2.8) 

(2.9) 

(2.10) 
(2.11) 



The following two lemmata, to be proved in §5.5 and §5.6, respectively, resume the main 
properties of -R($). 

2.8. Lemma. There are two constants C\ > 1 and Ci, independent of i, p and N, such 
that, if$,$'eB and 

dQi; - 4 |cr|[l + log( Wo 2 /|cr|)] <i, (2.12) 



then 



\\R(*)-R(&) 



< 



C 2 3 N Nl 



1 + log ■ 



v 



(2.13) 



2.9. Lemma. There is C > 1, such that, if 

CQ Uo -V| 1/2 [l + logK 2 /M)]<l, (2-14) 
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then 



P(0)||^<^fl + logl| sup 

«0 V v 0/ |n|>l 



(2.15) 



2.10. Lemma. There are e,c,K, independent of i, p and N, such that, if a £ J 

X 2 <r ^(l + logV)" 1 

" K N N\]og(cQ/i4) ' 
there exists a unique solution $ € B of (2.3); moreover the solution satisfies the bound 



(2.16) 



m-F< i — 



(2.17) 



2.11. Proof of Lemma 2.10. It is easy to see that, if a e J, the conditions on cr of Lemma 
2.8 and Lemma 2.9 arc satisfied, if 



A 2 <e /log(cQ/<) 



(2.18) 



with suitable values of eo and c. Moreover, if eo < £iVo(l + log-u^ 1 ) -1 an d £i is chosen 
small enough, (2.6) and (2.18) imply that A 2 |c„(ct)| < 1/2, so that, by using (2.7), (2.8) and 
Lemma 2.8, we have that, if $ £ B, 



|T A ($)||^<4A 2 (l + A 2 - 

vo 



ll^(0)||^- + — f 1 + log — ) 3^^111*11^ 

vo \ v 



(2.19) 



Therefore, by (2.15) and (2.4), there exist constants C3 and C4, such that, if e\ < 
svo(C^ TV!)" 1 and e is small enough, 



|T A ($)||^<^(l + logl 
v o V ^0 



iN AT , , 1 iJV ( 7rv o\n\\ 
3"Nl + sup |n| JV exp - 

!„!>! \ 1UA Z / 



< 1 . 



(2.20) 



Moreover, by (2.13), if and similar conditions on A are satisfied, we have 

||TA($)-TA($ ; )lk< 5 2 1 + log- ||<i>-$'||^< (2.21) 

V Q \ VqJ I 

The bounds (2.20) and (2.21) imply that B is invariant under the action of T A and that 
T A is a contraction on B. Hence, by the contraction mapping principle, there is a unique 
fixed point $ of T A in B, which can be obtained as the limit of the sequence <3?( fe ' defined 
through the recurrence equation <J>( fc+1 ) = T A ($( fe '), with any initial condition e B. If 
we choose = 0, we get, by (2.21), 



00 00 

|*||^ < £ - < £ ^zrll* (1) ll^ < ll<f (1) ll^ 



(2.22) 



i=l 



On the other hand, by (2.15), 



l<f (1) l^=l|T A (0)^<^^fl + logl) 

\ v J \v 



•» 2 



(2.23) 



which immediately implies the bound (2.17), if e\ < £Vq{C§ N\) 1 , with e small enough. ■ 

2.12. Let us now consider the equation (2.5). We want to prove that it has two solutions 
of the form (1.19), if cr G J and Li is large enough. In order to achieve this result, we need 
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some detailed properties of the function c\ (a) , which are described in the following Lemma 
2.13, to be proved in §4.10. We need also the bounds on pi(a, $(cr)) and its derivative with 
respect to a, contained in Lemma 2.14, to be proved in §5.7. 



2.13. Lemma. There is a constant C , such that, if 

v _ 1 
— ; — r < e < — , 

Li (7 ~ ~ 87T 



< 1 



then 



with 



-ci(ct) 



1 



log-^ +n{a) 



2ttv 

\ri(a)\ <C[ l + lo 



1 

"o 



dnja) 
do 



< C 



1 



(2.24) 



(2.25) 



(2.26) 



2.14. Lemma. If a G J, X satisfies the inequality (2.16), with e small enough, $(cr) is the 
solution of the equation (2.3) described in Lemma 2.10 and 



^ 2itvop 1 (<j,*(o)) 



then there is a constant C, such that 

\r 2 (a)\<C 
dr 2 (a) . C 



1/4 



do- 



< 



1/4 



(2.27) 



A2 

Vq 



(2.28) 



2.15. Lemma. There exist positive constants e, e, c andK, independent ofi, p and N , such 
that, if A satisfies the inequalities (1.16), there are two solutions a^(X) £ J of equation 
(2.5) of the form (1.19). 

2.16. Proof of Lemma 2.15. By using the definitions of r\{u) and r2{<j) given in (2.25) and 
(2.27), we can write the equation (2.5) in the form 



F{a) = log ■ 



+ r(cr) = 



(2.29) 



where r(er) = ri(er) + r 2 (cr). 
Let us now suppose that A satisfies the inequalities (2.16) and that a belongs to the interval 



J 



^2 e -4™ /A 2 ) v 2 e -nv o/A^ c j 



(2.30) 



If Li is large enough and the constant e in (2.16) is chosen small enough, the conditions 
(2.24) of Lemma 2.13 are satisfied, for a e J, and 



"A2" 



< log(ew Lj) 



(2.31) 



Moreover, if e and e (hence \cr\v 2 ) are small enough, 



dr(a) 



da 



1 

< - 
- 2 



d vl 

TT- log — 

oa a 



(2.32) 
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hence F(a) is a monotone decreasing function of a in J. If wc define 

a* = vle~ 2wvo / x2 , M = sup |r(<7)| , (2.33) 

we have that F(a* exp(-2M)) > and F (a* exp(2M)) < 0. Moreover, the interval 
(a* exp(— 2M)), c* exp(2M))) is contained in J, if e is small enough, since the bounds (2.26) 
and (2.28) imply that M < C(l + logw _1 ). Hence there is a unique solution <rW(A) of (2.29) 
in J, which can be written as 

cr (+) (A) = wge xa j (2.34) 

with \f3 {+) (\)\ ^CA^l + logi-- 1 ). 
In the same manner, we can show that there is solution cr(~\\) in the interval 

(-v 2 e-^ , -vle'^) C J , (2.35) 

with the same properties. ■ 

2.17. LEMMA. The constants e, e, c and K, appearing in (1.16), can be chosen so that the 
Hessian matrix (1-15) is positive definite. 

2.18. The proof of Lemma 2.17 is in §5.8. This completes the proof of Theorem 1.6. 
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3. Graph formalism 

3.1. In this section we shall describe the expansion of p x {ip,n), used to get the results of 
this paper. 

Let us consider the operators -0x = e tH ip x e~ m , with x = (a;, t), —(3/2 <t< (3/2 for some 
(3 > 0; on t antiperiodic boundary conditions are imposed. As explained, for example, in 
[BGM] , there is a simple (well known) relation between p x (ip, p) and the two-point Schwinger 
function, defined by 

qL 0i , Tr [exp(-0H) (9(x > y )^ x ?/>+ - 0(x a < y )V>x V^)] 
b ' (X;y) " Tr [exp(-/3ff )] ' 

given by 

p x = -\im lim \s LJ3 {x,t;x,Q) . (3.2) 

In [BGM], which we shall refer to for more details, it is also explained that the two-point 
Schwinger function can be written as 

^ ( x; y)=i im j ;;^ )eV ;f y , (3.3) 

where ^ are now anticommuting Grassmanian variables and P(di(j) is a Grassmanian Gaus- 
sian measure, formally defined by 

*>(#) = { II W(k))}exp|- (WO) -1 ^^ I #"# + , (3-4) 

k = (fc,fc„), 2? L ,^ = P L x Vp, V L = {k = 2irn/L,n G Z,-[L/2] <«<[(£- l)/2]}, 
D/3 = {fc = 2(n + 1/2)tt//3, n e Z, — M < n < M — 1}, in the limit M -> oo, 

— zfco + cospp — cos A; 

is the propagator or the covariance of the measure, pf = npis the Fermi momentum, defined 
so that cos p F = 1 — p,, and 

V(V) = V / <te A^V+C • (3-6) 

If we insert (1.10) in the r.h.s. of (3.6), we get 

I(Q-i)/2] 

v w= E ™ E A ^^c +2 „p F , (3-7) 

n=-[Q/2] M k€X>i,/s 

where pf = (pf, 0) and fc + 2«pf is of course defined modulo 2n. 

3.2. Note that g(k) _1 is small for k ~ ±pf- Hence there is no hope to treat perturbatively 
the terms with n = ±1 and k near =FPf, but we can at most expect that the interacting 
measure is a perturbation of the measure (whose covariance is not singular at k = ±Pf) 



P X (d^) EE ^P(#) 



eXP l A ^4 E \ E [^k+2p F +^k+2p^k] > 

[ M k av, 3 fee/_ J 



(3.8) 
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where AT is a normalization constant and 7_ is a small interval centered in — pp, so small 
that J_ n 1+ = 0, if 1+ = {k = k + 2p F , k e /_}. 

This remark suggests to apply a multiscale expansion to the integral (3.3), in order to 
treat in a different way the momenta near ±pp and the others. This procedure was applied 
in [BGM] to study systems of electrons in presence of a potential of the form (p(2px), with 
ip 27r-periodic, p/ir an irrational Diophantine number and pp = mp, m arbitrary integer. 
In [BGM] the aim was mainly to get the best possible results about the dependence of the 
two-point Schwingcr function on A and we found it useful to realize the multiscale expansion 
by dividing the momenta near pp into a number of "slices" of order | log A|. 

This expansion could be applied also the case pp = p with p = pp/tt rational, without no 
important difference, and we could get immediately Lemma 2.2. However, in this paper we 
prefer to use a simpler expansion into only two scales; this expansion gives weaker results 
about the dependence on A, but it is sufficient in order to prove Lemma 2.2 and it is more 
suitable for studying the equation (1.12). 

3.3. Let us introduce a smooth positive function fo(k') on the one dimensional torus T 1 , 
such that 

*<*>-{;: * K>^ 2, (3 - 9 ' 

where 

to =mm{p F /2,(TT-pp)/2} . (3.10) 
and the norm ||fc'|| T i on T 1 is defined so that H&'Hx 1 = W\, if k' € [—n,n). Then we write 

f 1 (k) = l- f Q (k+ PF )- fo(k-pp) , /o(fc) = l-A(fc), (3.11) 
so that (3.5) becomes 

m - »<■>« + = E z^S_ . ,3. 12) 



and, for h = 0, we define 
where, if k' = k — ujpf, 



(0) 



(k) = ]T 5i 0) (k) , (3.13) 



w=±l 



(k' + upf) = g £» (k0 = — — /0(fc/) ... 

— ik + cospp — cos(fc' + LOpp) 



(3.14) 



with k' = (fc',fc ); we set also ff^(k') = g (1) (k), with k = k' +p F , and /i(fc') = /i(fc), in 
order to simplify the notations in the following sections. 

3.4. We can associate with the decomposition (3.12) of g(k) a decomposition of the Grass- 
manian Gaussian measure P(dip) into a product of two independent Grassmanian Gaussian 
measures: 

P{dip) = P(d^)P{d^) , (3.15) 

if P(dV> w ) is defined as m (3.4), with #«(k) in place of g(k). 

If we insert (3.15) in (3.3) and we perform the integration over the field it is easy to 
show (see [BGM], §4) that 

5 L ^(x;y) = . 9 ( 1 )(x;y)+^(x;y)+5(°)(x;y) , (3.16) 
where g^fry) = (i/?)" 1 E ke ^,, 9 (1) (k) exp[-ik ■ (x - y)] and 



J dx(0jvi OH +V-i o)+ fc) V( o )(V <o) )+W (o) (^ (o) ,0) 



(3.17) 



<p+ = 
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in (3.17) J dx is a shortcut for J^xeA 

A/o = / P(dip( ^) exp[V(°''(V^ ' ) )], {0x} are Grassmanian variables anticommuting with {V'xl 
and V^ ** {fp^), the effective potential on the small momenta scale, can be easily represented 
as a series in A, as well as the function W^> 4>) and the function K^(x;y) appearing 
in (3.16). 

A precise description of these series in terms of Feynman graphs can be found in [BGM] ; see 
in particular the equations (4.3), (4.4) and (4.6) in [BGM]. Here we want only to stress that 
the involved graphs are chains formed by vertices, associated with the Fourier components 
of the potential ip n , connected through lines, associated to the propagators g^; hence by 
using the bound (3.27) below, it is easy to prove that these series are convergent, uniformly 
in L and (3. On the contrary, the series obtained by integrating the field ip^ do not have 
this property, for the reason explained before, and we have to look for a different expansion, 
based on the idea outlined in §3.2. 

We can associate with the decomposition (3.13) of g^(k) a decomposition ip^ = + 

V4°- °f t ne ne ^ ^k°^- The support properties of fo(k'), see (3.9), and the definition (3.14) 
imply that the field i/4°i nas support on the set {k = k' + Lupp ■ fo{k') ^= 0} and that the 
supports of ^4°+ and ^4°- are disjoint. 

The idea is to modify the free measure P(dip^) by multiplying it by the terms present in 
y( )(^j( )) ; (see [BGM], §3), which couple the variables v4°- and V4+2 PF +> tncn we expand 
the integral by using the new measure as the free measure. The new graphs differ from 
the previous ones for two respects; first of all they are not singular anymore at k = ±Pf, 
but they are bounded by C| A«^i | ~ 1 , see below; moreover the two propagators exiting and 
entering in the same vertex can not have both the momentum equal to ±Pf- As we shall 
see, these two properties are sufficient to control the expansions. 

The two properties of the new free measure described above are realized also if we only 
extract from V^(ip^) the first order terms coupling ^4°- an d V4+2p F +• ^ is easv to see 
that these terms are equal to \ip\F^ (ip^), with 

^ (0) )=EA E ^t PF ,.^_- PF; _. • 0.18) 

Hence we define 

P(^(o) } = ^ P(# (0) )e A^F(°>(v,<°>) ; (3 19) 

where J\f is a suitable constant, and 

y(0)ty,(0)) = v(°)(V (0) ) - A^Fi°)(^°)) . (3.20) 

By proceeding as in [BGM], §3, one can show that the Grassmanian integration P(dtp^) 
has propagator 

w,w'=±l 

with 

ffi°L'(x;y) = ^ E e - 4k '- (x - y) ffi°L(k') , (3.22) 



- ( o) m = H^Mi^W) 

9u ' uK ' l-iko-F^k'W-lFiW + aKk') 
M*')]/o(*') 



(3.23) 



L -ik - F^k')} 2 - [F^(k') + 4(k')\ ' 
CTo(fc') = crfo(k'), if a = X(fi ^ 0, and (see (1.17) for the definition of v ) 

Fi(fc') = (cosfc' - l)cosp F , F 2 {k') = v sinfc' . (3.24) 
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3.5. Remark. Note that, if \k'\ < t , 2(1 - cosfc')/| sinfc'| = 2|tan(fc'/2)| < 2tan(i /2) < 
| tan |; hence 

\Fi(tf)\ < \\F 2 (k')\ , for < t - (3.25) 
It immediately follows that 



ig(o) (k ' ){ < c V«o + W>«T fju) 

W k)| - 6 ^w 2 +/* ) ' (3 26) 

where C denotes a suitable constant. From now on, for simplifying the notation, the symbol 
C will be used everywhere to denote a generic constant, that we do not need to better 
specify. 
It is also easy to prove that 

|gW (k ')| < C - lzMt2 = , (3.27) 



3.6. We now insert (3.19) and (3.20) in (3.17) and represent the result of the integration 
in terms of Feynman graphs, by using P(dtp^) as the free measure and V(°)(^(°)) as the 
effective potential; then we apply (3.2). By proceeding as in [BGM], it is easy to show that 
we get an expansion for p x , which can be described in the following way. 

3.7. A graph $ of order q > 1 is a chain of q + 1 lines £\, . . . ,£ g +i connecting a set of q 
ordered points (vertices) v\, . . . ,v q , so that £j enters Uj and £ i+ i exits from Vi, i < q; the 
lines l\ and £ q +i are the external lines of the graph and both have a free extreme, while the 
others are the internal lines; we shall denote int(i?) the set of all internal lines. We say that 
Vi < vj if Vi precedes Vj and we denote v'j the vertex immediately following Vj , if j < q. We 
denote also by £ v the line entering the vertex v, so that £i = £ Vi , 1 < i < q. We say that a 
line £ emerges from a vertex v if £ either enters v (£ — £ v ) or exits from v (£ = £ v i). By a 
slight abuse of notation, if v = v q , we still denote by £ v > the line £ q+ \ exiting from v q even 
if there is no vertex v q+ \. 

We shall say that d is a labeled graph of order q > 1, if $ is a graph of order q, to which 
the following labels are associated: 

• a label n v for each vertex, 

• a frequency (or scale) label hi for each (internal or external) line, with the constraint that, 
if n v = ±1 for some v, then hi v — hi , = is not possible, 

• for each line £, two labels u>j, u>g, such that u>j = uJf = 1 if hi = 1, 

• a momentum — k = k' + uo\Pf for the first line, 

• a momentum 

ke v =k' + J2 [ 2n -PF + ( w l - w i )Pf] ( 3 - 28 ) 



for each internal line, 
• a momentum 



K+i = k ' + S [ 2n *PF + {"i ~ W/ 0/ )pf] (3.29) 



for the last line. 



If ij^ji u 2 (M-i) denotes the propagator associated with the line £, we will use the shorthand 
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Let us call T n>q the set of the labeled graphs of order q and such that 



2n v pp + (u)f — u>Y) pp = 2npp mod 2ir . 

veti lei) 



(3.30) 



Then, if p n (ip,p) is defined as in (1.11), we have 



Pn(V,M) = i im 



1 X 



L/3 



9=1 



p«(<7,4>) = X Val(tf), 



(3.31) 



where 

vaiw=~ e (n&)(n 

M k'er>i,, fl i=i veti 
Hence, the function p n (a,$) defined in (2.1) can be written as 



(3.32) 



p„(<7,$) = lim ^pl{a,$>) 



9=1 



#(*,$)=p«(a,$), if Q > 2 
^(a,$)= ]T (l-*n,n„)Val(tf). 



(3.33) 



after substituting in the r.h.s. of (3.32) A<p nt) cither with $ nt) , if |n„| > 1, or with ct, if 
K| = 1. 
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4. First order graphs 



4.1. In this section wc study the first order contributions to the density, i.e. the terms 
corresponding to graphs with only one vertex in the perturbative expansion (3.31), calculated 
in the limit (3 — > oo. For these graphs we have, if L = Li = iQ, 



lim Val(tf) 



A ^Z E r^^( k ')^(k' + (2m + ^-a; 2 )p F ) 



(4.1) 



where T>' L is the set of possible values of the variable k' introduced before (3.14) as the differ- 
ence between the "space momentum" k and ±Pf- Since pf = = ttP/Q = (2ir/L)(iP/2), 
we have 

V L = {k' = ^(n + 5/2), n G Z, -[L/2] < n < [(L - l)/2]} , (4.2) 

where <5 = 1, if iP (the number of particles) is odd, while 5 = 0, if is even. 

Note that the value of 5 will be in general not relevant, except in the proof of Lemma 2.13 
in §4.10, the only place were there is a non trivial dependence on the volume. 

Note also that, if the graph value (4.1) contributes to /3 ra (y, /it), then 



2mpp = 2npp + (wi — u>i + ^2 — w 2 )f>F m °d27r 



(4.3) 



The r.h.s of (4.1) can be easily bounded, by using (3.27) and (3.26) and the remark that 
lim PF ^o to/vo = 1/2. If h = h' = 1, one gets, for any integer r, 



fc'GP'r 



C / dfc' 



dko 
~27 

dko 



g[%')g[%' + 2rp F ) < 
[l-/o(fc')] [1 - /o(A;' + 2rpF)] 



+ («ofc') 2 v^o + W + 2 ^) 2 
If h = 0, /i' = 1, for any w, u/, one gets 

i e r$|e(k')^(k' + 2, P ,)| 

^ h c tv •'-oo 27r 



< 



(4.4) 



< 



C 



o 2tt y v /fc^ + ( Uo fc') 2 yfcgT^ 



c_ 

v 



(4.5) 



The bound (4.5) can be improved for a — > 0, if a/ = — w; we have 



dk 
2tt 



fc'ex>' r 



~9 i "l OJ (k')~g[':i(k' + 2rpF) 
dkp 



< 



[k 2 + (v k') 2 + c7 2 }^/kj+^ ~ v$ 



< 



1+log n) 



(4.6) 



Let us now consider the case h = h' = 0; for any ua, w-, we get 



L t-cW 7 -o« ^ 

/o(fc') 



fc'eiV 



< 



C 



C 



dfc' 



dko 

> 

/■to /■oo 
JO JO K 



f (k' + 2rp F ) 



^k 2 + {v k') 2 + a 2 V^o + v o( k ' + 2r PF) 2 + <? 2 

* ^<^i + i» g ^). 



< 



(4.7) 



+ {v k') 2 +, 
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If h = h' = and ui\ ^ u[ or uj 2 ^ u)' 2 , the last bound can be improved; in fact we get 



< 



k'EV' 



C 



a\ / dk' 
Jo Jo 



dko 



< 



C 



(4.6 



[k 2 + (v k') 2 + a 2 }^/ 2 " «o ' 



The previous bound can be further improved, if we suppose also that r ^ 0, by taking 
into account that, in this case, max{|fc'|, \k! + 2rp F \} > n/Q. Let us suppose, for example, 
that 1V2 = —u>2, we have 



- e r 



dko 
~2tt 



dk' / dk a 

-IT JO 

c\*\q r dk , ■ 



fo(k') 



< 



fo(k' + 2r PF ) 



^k 2 + (v k') 2 + ° 2 [k 2 + v 2 (k' + 2rp F ) 2 + a 2 } 



< 



dko 



C\a\Q 

vo Jo """ Jo k o + ( v ok') 2 + <t 2 ~ v 2 



< 



1 + log ft 



(4.9) 



In the following we shall need also the bounds of the expression obtained substituting in 
the r.h.s. of (4.1) one of the two propagators with its derivative with respect to a. By 
proceeding as before, one can easily prove that, for any to and any integer r, 



L 



k'GT>' T 

dko 
~2tT 



dko 
~2tt 



dgffL(k') _ ( i) , 



da 



5i^(k' + 2rp F ) 



< 



C_ 

IS > 



d£°L(k') _(i) , 



9a 



<§(l + log"' 



that, for any and any integer r, 

and finally that, for any u>\, u[, oj and any integer r^O 



< 



c 



■ E 

fc'eix 



00 &o 

-oo 27T 



<9er 



£ 0) _Jk' + 2rp F ) 



j. y°° dko 

L y^oo 2tt 



,% (0) 

CU( k ')^( k '+ 2 ^) 



< 



< 



CQ 

2 ' 



CQ 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



4.2. Remark. All the previous bounds are valid also if we exchange in the l.h.s. k' with 
k' + 2rp F ; this immediately follows from the observation that the variable k' is defined 
modulo 2n. 

4.3. We shall now consider the graphs contributing to the constants c n (a) introduced in 
(2.1), in order to prove Lemma 2.6. We can write 



-\(p n c n {a)= ^2 <W„Val(?9) • 



(4.15) 
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The equations (4.15), (4.1) and (4.3) imply, if |n| > 2, 

c»4 E r^{^(k')^(k' + 2np F) 

+ E (k') 5i°L(k' + 2np F + (1 - w )p F ) (4.16) 

«=±i 

+gi°l(k').g«(k' + 2np F - (1 - w )p F ) + £°L(k') .g^(k' + 2np F ) 
+£°L(k') 5i°L(k' + 2n PF ) + gW(V) £°>,_(k' + (2n + 2u;)p F )] } . 

By using the bounds (4.4), (4.5) and (4.8), we see that the first four terms in the r.h.s. of 

(4.16) are bounded by (C/t>o)(l + logw^" 1 ). However, the remaining terms, i.e. those with 
h = h' = and u)\ — ui[ = uj 2 — u>' 2 = 0, need a more careful analysis; these terms will be 
denoted as 

°» s 7 E T ^£°L(k')£°L(k' + 2n PF ) , (4.17) 

when loi = uj 2 , and 

b n ,„ = }E r ^ ^l(k') £°2,_ w (k' + (2n + 2u;)p F ) . (4.18) 

when uii = —u 2 . The following two Lemmata 4.5 and 4.7 show that the dimensional bounds 
which would follow from (4.7) in fact can be improved. 

4.4. Remark. Note that a n is a ^-independent quantity, so that we can set ui = 1 in 

(4.17) ; this property easily follows from the observation that g£°L(fc', k a ) = g^ t _ u (—k' , fc ), 
see (3.23). It is also easy to prove that b ni i = b- n> -\. 

4.5. Lemma. Let \n\ > 2 and let a n be defined as in (4-17); then \a n \ < C/v . 

4.6. Proof of Lemma 4-5. By Remark 4.4, it is enough to study the case ui = 1 in (4.17). 
Define 

=(0) / k n = fo( k ') 

■ ] -iko ' (4.19) 



and 



T(k') = sign(k')^Fi(k') + a 2 (k') - F^k') , 
MjE r^(k')^(k' + 2np F ) 



(4.20) 

= L E /o(fc').fo(fc' + 2n^)X(fc') • 
k'ev' L 

where A n (k') is obtained by explicitly performing the integral on ko- It is easy to see that, 
defining 

s(k') = sign (>(fc')) , (4.21) 
if s(fc') = s(fc' + 2np F ), one has ^4„(/c') = 0, while, if s(k') = —s(k' + 2np F ), one has 

An(k') = s(k') \j=-(k' + 2np F ) - F{k')\ 1 . (4.22) 

Note that, by (3.25), s(k') — sign(fc'), if \k'\ < t , i.e. on the support of fo(k'); hence we 
have 

- _ J_ h{k')f (k' + 2n PF ) 

n ~ l k , e ^ nv , + \Hk' + 2n PF )\ ' [4 ' 26) 
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where 

£>* = {*' G [-to, h] ■ sign (k 1 ) = -sign (k 1 + 2np F )} . (4.24) 
We want to show that 

max{\F(k>)\, \nk' + 2n PF )\} > ^-\\2np F \\ Tl = Ai , (4.25) 

if k! e T>[, k' + 2np F e [-t ,t ] and c 2 = (V2/ir)v . 

If \T{k')\ > Ai, (4.25) is immediately verified. Let us suppose now that |-F(fc')| < Ai; 
then, by using (3.25), we get 

\Hk')\ > \F 2 (k')\ \ Fl (k')\ > \\F 2 (k>)\ > c 2 \k'\ , (4.26) 

so that 

< = l\\2np F \\ Tl , (4.27) 

C 2 Z 

implying 

\\k' + 2np F \\ Tl > |||2np F || Tl > ^||2np F || Tl . (4.28) 

Moreover, since \\k' + 2np F \\j 1 < t n , then 

\T{k' + 2n PF )\ > c 2 \\k' + 2np F || Tl ; (4.29) 

hence, by using (4.28) and (4.29), we get 

\T{k' + 2np F )\ > ^\\2n PF \\ Ji = A 1 , (4.30) 

which implies (4.25) also when \T{k')\ < Ai. 
Inserting (4.25) into (4.23) leads to 

K < / ^ jp- < — , 4.31 

J v , r nc 2 \\2np F || Xl v 

as the size of the set T>[ is bounded by 2||2npF||Ti • 
In order to complete the proof of Lemma 4.5, we note that 



J E r^{[^(k')-^(k')]^(k' + 2npF) 
+ g[°l (k') [g[°l (k' + 2np F ) g[°l (k' + 2np F )] } . 



Moreover, by (3.23) and (4.19), 



(4.32) 



|-ifc - Fi(fc')] 2 - [F4{k') + a${k')]\ 



a„-a n \<C dk' / rffcp - < — . (4.34) 

Jo io v/fwofc') 2 + o 2 \kl + (v k') 2 + a 2 ! 3 / 2 " v 



so that, by using also (3.25), we get 

fto f-oo 2 

^(Wofc') 2 + [jfcg + (W fc') 2 + ^ 2 ] 3/2 ~ W 

The bounds (4.31) and (4.34) imply Lemma 4.5. ■ 

4.7. Lemma. Let \n\ > 2, and let b n ^ be defined as in (4-18); then \b n ^\ < (ClogQ)/i> . 

4.8. Proof of Lemma 4-. 7. Let us define b n ^ as 



b»„ = \ E r tflWS^-vKM , (4-35) 
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where k' u = k! + (2n + 2oj)pp and 5^°L(k') is denned in (4.19), and set b ntU! = 6„ ;W + J'. 
By dimensional bounds analogous to those which led to (4.34), it is easy to prove that 
| J' | < C/vq. Moreover, by proceeding as in §4.6, we see that 

h -I ^ /o(t')/oK) 

where = {k' e [-io,*o] = sign(wfc') = sign (wk£,)}. 
By using the bound (4.26), we have 

|^fc')l + > cadfc'l + IKHtJ (4-37) 

Moreover, since |n| > 2, ||2npF ± 2pf||ti — ~£?> nence 

- /■ dk' ClogQ 

- 27t C2 (i*'|+ikii Ti ) - • (438) 

This completes the proof of Lemma 4.7. ■ 

4.9. Proof of Lemma 2.6. The bound (2.6) immediately follows from the remark after 
(4.16), Lemma 4.5 and Lemma 4.7. The bound (2.8) is easily proven from (4.16) by using 
the bounds (4. 10)-=- (4. 12). ■ 

4.10. Proof of Lemma 2.13. The definition of ci(<r) in §2.3, (3.31), (4.1) and (4.3) imply 
that 

l ^ r°° dk f~g%(V) (1) (1) 

k'ev' L 00 

+ ^ [.g«(k0.g£%(k' + (3- w )p F )+.9^(kO^(k' + (l + 

w=±l 

s-F^.LJ + ci^), (4.39) 

where —F(a,L) denotes the first term in the r.h.s. of (4.39), while c\(a) is the sum of the 
other ones. It turns out that — F(<t, L) is the leading term for a — > 0; moreover it is the only 
term whose dependence on L is not trivial, hence we decide to indicate it explicitly. 
By using the definition (3.23) and by performing the integration over fc , we get 

F(a, L) = J- £ - /o(fc,)2 . (4.40) 

2L k>ev' L ^ 2 sin 2 fc' + ( 7 2 /o(fc') 2 

The definition (4.2) of T>' L implies that, for any finite volume L = Li, the r.h.s. is singular 
for a — > 0, only if S = 0, that is only if the number of Fermions is even; in that case, in fact, 
k' = belongs to the set T>' L . It follows that, if 5 = 1, the equation (2.5) has no solution for 
A 2 very small, how small depending on L; this is the main source of the lower bound on A of 
Theorem 1.6. [Equivalently, for fixed A verifying the inequality to the right in (1.16), which 
is uniform in L, L has to be large enough so that also the inequality to the left in (1.16) can 
be fulfilled.] We separate the term with k' — 0, if it is present, by writing 

F(a,L) = ^-+F (a,L) . (4.41) 

It is easy to see that 

F {a,L) = F 1 (a,L)+d 1 {a,L)+d 2 {a,L) , (4.42) 
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where 



FML) ^o^-IWv* V(2-o) 2 (n + S/2) 2 + (aL) 2 



A2 



k>zv' L yjv%sm'k' + a 2 f (k') 2 

\k'\>t /2 



0^|fc'|<t /2 



W 2 sin 2 jfe' + a 2 V(vok') 2 +<r 2 



(4.43) 
(4.44) 

(4.45) 



Note that the sum in the r.h.s. of (4.43) is empty, if [toL/(4n) + S/2] < 1; in that case 
the equation (2.5) may have a solution, for A small enough, only if S = 0. Hence we shall 
suppose that: 

t L > 4tt , (4.46) 



a condition which is certainly verified, if the conditions (2.24) are satisfied, since 

- < v /t < 2 . 



(4.47) 



By using (4.47) and supposing that 



it is easy to show that 



Ei*(*)i^. E 

i=l i=l 



ddi(a) 



da 



< 



C 



By substituting the sum in the r.h.s. of (4.43) with an integral, we can write 

F 1 {a,L) = F 2 {a,L) + d 3 (a,L) , 

where 



F 2 ( 



t L/(4n) 



dx 



-8/2 y/(2TTV x) 2 + {<TL) 2 

It is easy to see that, if the condition (4.48) is verified, together with 



then 



< s < 1 



dd 3 (a) 



da 



< 



Ce 
vo\a\ 



The integral defining F 2 (a,L) can be explicitly calculated; we get 



F ^ L) =2^ l ° g - 



vgta 
2\a\ 



S(i-i)+Vte)"( 1 -^) a+1 



If we write 



F 2 (a,L) = -^— log ^L + d 4 (a,L) , 
2irv \a\ 



(4.48) 

(4.49) 

(4.50) 
(4.51) 

(4.52) 
(4.53) 



(4.54) 



(4.55) 
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it is easy to prove, using (4.47), (4.48) and (4.52), that 



ddi{cj) 



da 



v 



(4.56) 



Finally, the function Ci(cr) introduced in (4.39) and its derivative can be bounded, by 
using (4.4), (4.5) and (4.10), as 



ci{a)\ < — ( 1 + log — 



It is now sufficient to define 

r\(a) = 2nvo 
to complete the proof of Lemma 2.13. 



vo 



dci (ex) 



da 



C 

^ — ■ 

v o 



1-6 

Y J d i {a,L) + ^—-c 1 (a) 



(4.57) 



(4.58) 
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5. Bounds on the density perturbative expansion 



5.1. In this section we give some bounds about the perturbative expansion (3.33) of the 
function p n (cr, $), introduced in (2.1), and we prove Lemmata 2.2, 2.8, 2.9 and 2.16. 
Given $ £ T , let us define 



R($)M(a)=pUa,$(v)) , M>0, q>0. 



(5-1) 



Moreover, if J is the space of the C 1 -functions of a £ J with values in R, and r(er) £ J, 
we shall define, in agreement with (2.9), 



11711.7 = sup 



W\- 1 \r(*)\ + 



dr 



(5.2) 



5.2. Lemma. // $ £ B and \a\ < Vq, then, for any n^O and q > 0, 



\\R{*)% 



(9)1 



J 



< 



c 



3q 



M 1+log WUJ q \™ 

where C and D are suitable constants. 



N 



1 + (1 log ^ 
\o~\ 



W\Q) 



[9/2] 



(5.3) 



5.3. Proof of Lemma 5.2. In order to bound p\(cr, we shall use the expansion in (3.33). 
Let i3 £ T n:q be one of the graphs contributing to p^(cr, $) and v one of its vertices. If 
\n v \ 7^ 1, one has (see §3 for notations) 



> 



2tt 

~Q 



so that 



max-tll^jT 1 . II^JIt 1 } > q ■ 
Then there is a constant C2 such that, Vw € 1?, |n„| ^ 1, if |<r| < 1, 

C 2 Q 



< 



min{ |ff/J,|5/„,|} < -V 



(5.4) 
(5.5) 

(5.6) 
(5.7) 



by (5.5), (3.27) and (3.26), since v < 1. 

Note that (5.6) and (5.7) still hold for |n„| = 1, as, in such a case, hi v — hi , — is not 
allowed (see §3.7) and the support properties imply that both propagators are bounded by 

C/vl 

Note also that, thanks to (3.30), 



M <g+l + ^K| <3 



3v* 



> 



3q 



(5i 



Let us now suppose that q = 2q, with q > 1. It is easy to see that, in this case, it is 
possible to couple 2q among the 2q + 1 propagators appearing in the expression of Val(#), 
see (3.32), in q pairs {ge v ,ge ,} with v ^ v*; let g#i), £^ = £„*, the propagator left alone 
after this coupling operation. We select in an arbitrary way one of the q couples and we use 
the bound (5.7) for one of the propagators belonging to it; let g t ( 2 ) the other propagator of 
the selected couple. The propagators of all the other couples will be bounded by (5.6). We 
get 

-9+1 . 2 « 



|Val(0)| < (C 2 Q) 



4q-2 



(n 

i=l 



(5.9) 
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Let us now suppose that \a\ < u 2 ; then we can use the bounds (4.4)-(4.8), valid also for 
finite (3, to prove that, for any choice of g^i) and g £ (2) , 



- E 



|&(1)<?£(2)| < — ( 1 +10g — 



1 



(5.10) 



Hence, if $ e B, by using (2.10), (5.8) and (5.10), we get 
E |ValWI<|^(l + lo g l) 

D?gCfM(MQ)«"fl + log^ 



#6T n ,2,- 



JV 



(5.11) 



2 



where D'g takes into account the fact that there are 5 possible choices for the uj,0J e , 
labels for each line, and q possible choice for the vertex v*; then the bound (5.3) is proved 
for even q. 

The case q = 2q + 1 , with q > 1 , can be treated in a similar way. We note that it is always 
possible to couple 2q among the 2q + 2 propagators appearing in the expression of Val(#) 
in q pairs {ge„,ge ,} with v ^ v*; let g tm and g^ 2 ) be the propagators left alone after this 
coupling operation. Then we use (5.6) for all the couples and the bound (5.10) for the two 
remaining propagators. We get a bound similar to (5.9), with | cr| ~ "? in place of |cr|~^ +1 , but 
the final bound is the same as before. 

We still have to consider the case q — 1 . We could of course get again the previous bound 
with q = 0, but there is now an improvement, which will play an important role. The 
improvement follows from the observation that, if q = 1, the graphs contributing to p\ have 
only one vertex with Fourier index n Vl =/= n, so that at least one of the two propagators must 
have different w-indices. By using the bound (4.8), this implies that the bound (5.10) can 
be improved by erasing the factor [1 + log(v 2 /|cr|)]. 

In order to complete the proof of (5.3), we have to bound also 
9p^(cr, $(a))/ da. We can proceed as before, by noticing that dVal(6)/da can be written as 
the sum of 2q + 1 terms, each term differing from Val(0) only because there is the derivative 
acting on a single propagator or a single vertex function. If the derivative acts on one of the 
coupled propagators, one can use the bounds (5.6) and (5.7) modified so that the r.h.s. is 
multiplied by | cr | 1 ; if the derivative acts on a vertex function, since <f> e B, one can use the 
bound \d<& n (<j)/da\ < \n\~ N ; if the derivative acts on one of the propagators left alone after 
the coupling operation, one can use the bound, following from (4.10)-j-(4.12), if \a\ < Vq, 



- E 



dg e 



1) - 
9ew 



We get, for any q > 0, 



(</) 



da 



■(a) 



< 



E 



dVal(tf) 



< 



va\a\ 



(5.12) 



da 



< 



(2q+l)^ T DfqCl(\a\Qr(^- 



N 



(5.13) 



with q = [q/2]. This complete the proof of (5.3). ■ 

5.4. Proof of Lemma 2.2. The bound (5.3) immediately implies that p^(cr, $) is summable 
over q, for aQ/vg small enough, uniformly in i, p and /?. On the other end, it is easy to 
see that the bound is valid also if we substitute in the expression (3.32) of Val($) the sum 
over k with the integral on the real axis and that lim^oo J2^Li P%k G i ^) i s obtained from 
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Yl^Li PrSFi *) by doing this substitution. The claim of the Lemma about the continuous 
dependence on A of p n (f, p) i s an eas Y consequence of this remark and (3.31). 

In a similar way, one can see that lim^oo p n (<p, p) is obtained by substituting in the 
expression of Val(z9) the sum over k' with an integral over the interval [—n, n] and that this 
limit is also continuous in A near 0. 

The other claims of the Lemma about the density and the gap of h can be proved as in 
[BGM], §4.5 and §4.6. In [BGM] a more complicated expansion was used (involving a further 
decomposition of the field ip^), but the proof of these two points can be even more simply 
reached by using the expansion of this paper and bounds of the graphs similar to (5.11). 
We shall not give here the details, but we only remark that the main point in the proof is 
the remark that the propagators (3.23) are analytic in ko in the strip |Imfco| < H/2. ■ 

5.5. Proof of Lemma 2.8. By the remark in §5.4, lim^oo can be exchanged with the sum 
over q in (3.31) and (3.33). In the following, for simplicity, we shall use the notation p1(a, $) 
to identify lim^oo p\ (a, $). 
Let us suppose that $,$'e8 and q>2; then, by (3.31) 

pU<r, ]T Val(tf), (5.14) 

#£r„,,(4>,$') 

where T n q ($,$') is a set of labeled graphs whose definition differs from the definition of 
T n q , see §3.7, only because there is a new label a v € {0, 1,2} for each vertex; moreover 

vai W =4 e /"iKn&Xn^-)' (5 - 15) 

!<r if a v = , 

$ nv ifa„ = l, (5.16) 

$^-$„„ ifa„ = 2, 

and there is the constraint that at least one vertex has label a v — 2. 

By proceeding as in the proof of Lemma 5.2 and using the definitions (2.9) and (5.1), we 
get, if q > 2, 

iii?($') (9) - R(^) iq) y < 

Dv (l + log 1) ||$ ' - $||^ q \3q) N (l + log j|) (lalQ)^ . (5 ' 17) 

Hence, if Q^ 4 |cr|[l + log(ug/|cr|)] < 1 and |cr|0^ 4 < 1/(2C 2 ), we have 

V ||i?($') (9) - R(<S>) {q) \\r < — (l + log — ) 3 N N\ ||$' - $||^ . (5.18) 

with a suitable constant C\. 

In order to complete the proof of the lemma, we have to estimate ||i?($') ( ^ 1 ^ — i?($)^^ 
The bound (5.17), with q = 1, is still valid, but it is not sufficient; however there is the 
improvement with respect to (5.17) due to the fact that, if i? is a graph contributing to 
/5* (<r, $') — (c, $), the only vertex belonging to d has a Fourier index n Vl ^ n. As in the 
proof of Lemma 5.2, this remark allows to eliminate the factor [1 + log(vo/M)] in the bound 
(5.17) for any value of n. The previous remark implies that 

p($')« - R($)W\\r < £ (l + log ^) II*' - *lk ■ (5-19) 
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This bound and (5.18) immediately imply Lemma 2.8. ■ 

5.6. Proof of Lemma 2.9. The graph expansion of p^(a,0) has the property that, given a 
graph $ 6 T n q with Val(i?) ^ 0, each vertex of •& has Fourier index n v = ±1. This implies 
that, for any net? (see §3.7), he v = hg v , = is not allowed, so that the number of non 
diagonal propagators is less or equal of q + 1, if q = [q/2]. Hence (3.30) implies that 

\n\<q + q + l. (5.20) 

We can bound Val(^) ^ as in §5.3, by choosing in an arbitrary way the vertex v* (since 
we do not need now to extract the factor |n| _Ar ), and we get 

E iValWI^^fl + logl^S^^laKHQ^fl + logfi) . (5.21) 
It is easy to see that, if q > 2, q > |n|/5; hence, if Qu ( 7 3 |o'| 1 / 2 [l + log(7jQ/|cr|)] is small enough, 

In order to complete the proof of Lemma 2.9, we have to improve the bound (5.21) in the 
case q = 1. Note that p^(cr, 0) is different from only if |n| = 2 (only one propagator may 
have frequency label h = 0) and it is given, if n = 2 (the case n — — 2 is similar), by 



I £ r ^^^-(O) (k > + 4pF)+ £0) (k ^(i) (k ' + 2pF)] . (5 . 23) 
Hence, by using (4.6) and (4.11), we get 



2/10 

ll^(0)2 i; IU<^H(l + logf^) <- l^f) • (5.24) 



\(i)m - C 











(9) 




«0 





This bound and the bound (5.22) imply (2.15). ■ 

5.7. Proof of Lemma 2.14- By Lemma 5.2, if Q7J ' 3 |cr| 1 / 2 [l + log(uo/|cr|)] is small enough, 
which is certainly true if condition (2.16) is satisfied, with s small enough, we have 

El ^, $)l , H ^( 1 + log l)(H)" a . ( , 25) 

^2 v V v J V v J 

Moreover, since the graphs contributing to pj(cr, $) have only one vertex with index n v — 
±2, ±3, at least one of its two propagators has frequency index h = and different w-indices. 
It follows, by (4.6), (4.8) and Lemma 2.10, that 

\p{(a,$)\<-(\$ 2 (a)\ + \$ 3 (a)\)<%\(-) . (5.26) 
vo v \v J 

Hence, if \a\ 1/4 C N N\(l + log^ 1 ) < v]/ 2 , which is certainly true if condition (2.16) is 
satisfied, with e small enough, and r 2 ((r) is defined as in (2.27), we get 



\~p\{o)\ <&- 

V 



which implies the bound in the first line of (2.28). 



1/4 

10-1 1 



(5.27) 
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Let us now consider the derivative of r 2 (c). By Lemma 5.2, if |cr| 1 / 2 [1 + log(wo/|cr|)] Qv 3 
is small enough, we have 



E 

q>2 



da 



< 



C N N\ 



1 + log ■ 



v J \vfi 



1/2 



(5.28) 



Moreover, by Lemma 2.10 and the remark preceding (5.26), 



dp\{aM<y)) 



do 



C / x 2x N 



< -\\*y< 



(5.29) 



It follows that, if |<t| 1 / 4 C jv ./V!(1 + logt> *) < v^ 2 , which is certainly true if condition (2.16) 
is satisfied, with e small enough, 



dr 2 
~da~ 



< 



2ttvq 



WW 



1\\J 



< 



c 



1/4 



(5.30) 



which immediately implies the bound in the second line of (2.28). 



5.8. Proof of Lemma 2.17. The Hessian matrix M, defined in (1.15), is a real matrix; hence, 
we have to show that 

I(Q-i)/2] 

^2 x n M nm x m > , (5.31) 

n,m=-[Q/2] 



for any { x n}n2~[Q/2] e 3 - This wm ^ e done by writing 

I(Q-i)/2] 

n,m=-[Q/2] 



[(Q-l)/2] 



E 



n=-[Q/2] 



M„ 



, I(Q-i)/2] 

£ ( 

m=-[Q/2] 



M„ 



|M„ 



(5.32) 



and showing that the right hand side of the above equation is strictly positive. 
Let us find first a lower bound for M nn . If |n| ^ 1, by (1.15), (1.12) and (2.1) we have 



M„ 



l + X 2 c n (o)-X 2 9pn 



(5.33) 



where 1 + A 2 c„(er) > 1/2, see §2.11, and dp n /d<& n obeys to the same bound of dp n /do, see 
§5.3, up to the factor |n| _JV : simply note that the derivatives can act only on the vertex 
functions (and not on the propagators), and \d$ n /do\ < \n\~ N has to be replaced with 
\d$ n /d<& n \ < 1. Then, analogously to (5.13), we obtain, for any q > 0, 



2d~Pl 



< X 2 Dv 



q 2 {M) N {\^\Q) [q/2] 



(5.34) 



so that, if |cr|<5w 4 is small enough, we have 



,2 d Pn 



< CX 2 Vq 1 3 N N\ . 



It follows that 



H * 1 , 



(5.35) 
(5.36) 
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for A satisfying (1.16), with e small enough and K > 3. 
In the case n = 1 (the case n — — 1 is discussed in the same way) we have 



Mn = 1 + A 2 ci(cr) + Act 



da(a) _ ^dpi 



(5.37) 



dpi dipi 

Note that our definitions of Ci(a) and /5„(<7, <&) do not distinguish the dependence on (pi and 
if>-i, which are equal in the fixed points we are studying (see discussion in §1.4). However, 
in the definition of M, ipi and <p_i have to be treated as independent variables. By taking 
into account this remark and by using Lemmata 2.13 and 2.14, with e and |cr|z;^" 2 small 
enough, we get 

0ci(<r) 1,2 da(*) ^ 1 A 2 
= ^A a— > , 

2 OCT 07T Vn 

cx 2 



Xa- 



d0i 



1 + A 2 ci(a) 



dpi 



X 2 ~ P l 



< 



9<pi 



< 



> 9pi 
da 



< 



v 

CX 2 



v 



1/4 



1/4 

'-i I +1 " 



A 2 
vo 

A 2 



N' 



(5.38) 



AT 



so that 



, 1 A 2 

M nn \ > - 

8tt v 



±1 



under the hypotheses of Theorem 1.6. 
The non diagonal terms (n ^ m) are of the form 

9c m (a) 5 n s + 



M n 



A 2 $ r , 



A 



(5.39) 



(5.40) 



da 2 d(p n 

By using (2.7) and (2.17), the first term in the r.h.s. of (5.40), where m ^ n implies \m\ > 1, 
can be bounded as 

dc m (a) £>„,! + 8 n -i 



A 2 $, 



da 2 

Moreover, by proceeding again as in §5.3, we obtain 



< C 



N+l 



dpi 



difin 



< X 2 Dv [ 4 



p\Q) 



[9/2] 



so that 



[(Q-l)/2] 

E E 

n=-[Q /2];njtm 9=2 



dpi 



d(p r , 



< CX 2 3 N N\ 



HQ 2 



(5.41) 



(5.42) 



(5.43) 



The contributions with q = 1 need an improved bound. Let us first suppose that \n\ > 1; 
in this case the derivative can act only on the vertex function of the graphs contributing to 
p] n . Then, if the derivative is different from 0, the vertex function is equal to $ n and, since 
to — n ^ 0, at least one of the two propagators must have different u indices; this follows 
from (4.3), which also implies that the integer which multiplies pf in the value (4.1) of the 
graph is different from 0. Hence, by using (4.6), (4.9) and the fact that \n — m\ < 2, we get 
KQ-i)/2] Q _! . . 



E 

-[Q/2];n^m 



■9 Pin 



< cx 2 



logf^ 



(5.44) 



The case q = 1 and \n\ = 1 can be treated in a similar way; the main difference is that the 
derivative can act also on the propagators of the graphs contributing to p^, but it is still 
true that the integer which multiplies pf in the value (4.1) of each graph is different from 0, 
an essential point in the previous bound, since it allowed to use the improved bound (4.9) in 
place of (4.8). By using again (4.6) and (4.9), as well as the improved bounds (with respect 
to (4.12)) (4.13) and (4.14), we get again the bound (5.44). 

The r.h.s. of (5.41), (5.43) and (5.44) can be made arbitrarily small with respect to X 2 /vq, 
by suitably choosing the constants in (1.16); hence Lemma 2.17 is proved. ■ 
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